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1 Introduction 



The inflationary Universe model proposed by Guth [1] and its improved version by Linde 
[2], Albrecht and Steinhardt [3], is based on the theory of phase transition in the grand 
unified theories [4]. According to Guth's scenario,, as the Universe expands and cools down, 
a phase transition occurs from the symmetric state to the asymmetric state. During the 
phase transition the energy of the Universe will be dominated by the vacuum energy of the 
Higgs field, and the Universe is in a de Sitter state. It is this exponential expansion that 
the Universe has with a scale factor of at least 28 orders of magnitude, which resolves many 
long-standing cosmological puzzles such as the horizon and flatness problems. However, as 
was noted by Guth himself, the original inflationary scenario would result in an extremely 
large inhomogeneous and anisotropic Universe. 

This problem can be avoided in the improved version [2, 3] in which the first-order phase 
transition in Guth's model is replaced by the Coleman- Weinberg [5] symmetry breaking. 
Because the Coleman- Weinberg potential is very flat near the origin, the Higgs field will roll 
down the effective potential very slowly and the Universe keeps expanding exponentially. It is 
during. this period of slow growth that a single bubble inflates sufficiently to cover the entire 
Universe which will then result in a very homogeneous and isotropic Universe as we observe 
now. The success of the improved inflationary Universe model has led many authors to study 
quantum fields in curved space [ 6] to see the gravitational effect on the cosmological phase 
transition [7, 8]. From the effective potential they calculate, it was found that the scalar 
gravitational coupling ( and the magnitude of the scalar curvature R crucially determine the 
fate of symmetry. One knows that, in the inflationary scenario, the isotropy of the Universe 
before the phase transition is an a priori assumption. Although the quantum field effects 
on the curved space could account for this, such as the isotropization by particle production 
[9], it is less clear whether an initially large anisotropy can indeed be smoothed out at a 
sufficiently early stage, Therefore it is interesting to investigate the anisotropic effect [10-12] 
on the Coleman- Weinberg symmetry breaking. For the same reason, it is also worthwhile 
to investigate the effects of inhomogeneity [13! and rotation [ 14] on the cosmological phase 
transition. 

In this paper we will analyze the process of Coleman- Weinberg symmetry breaking in 
a rotational spacetime of Godel metric [15] at finite temperature. The zero-temperature 
theory in Godel spacetime has been studied in [14]. Note that the finite-temperature field 
theory in flat spacetime was investigated several tears ago [16]. It has also been studied 
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in the static [17] and dynamic spacetimes [18]. (A more complete list of references on the 
finite-temperature field theory in curved spacetimes can be found in [18].) 

We know that the problem of obtaining a solution to Einstein's equation for the rotating 
spacetime is fairly old. Godel [15] was the first to find an exact rotating solution. However, 
this solution suffers numerous pathologies, such as the violation of causality, the absence 
of a global Cauchy surface, the non-orientability of time, etc (see the discussion in [19] for 
details). This solution also has no expansion and no shear. Although the Godel model is 
not very physical, we study it here because it appears as the simplest rotating cosmological 
model. Also, despite the fact that the cosmological rotation is very small by the present 
observations it cannot be completely ruled out, at least in the early Universe. 

The paper is organized as follows. We begin in section 2 by solving the scale wave equation 
in the Godel spacetime and setting a formula to obtain the effective potential which, however, 
is formally divergent. In section 3 we use the ( function regularization method [20] to render 
it finite. A short discussion and conclusion are given in the last section. 

2 Formalism 

The action describing a massive (m), self-interacting (A) scalar field (0) coupled arbitrarily 
(xi) to the gravitational background is 

S[$] = J d'x^ {±<Td,M>4> - ^ - l R P ~ ^) ■ 2.1) 
The effective potential V e ff is given by [6] 




where A = g^V u and M 2 = m 2 + £R+ |A0 2 , with =< >. The above equation can 
be rewritten as 

J d 4 x^/=g~V eff = -^iTrln(A + M 2 ). (2.3) 
The trace can be easily evaluated once we know the eigenvalue of the operator (A + M 2 ). 
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Denoting the eigenvalue by r] and eigenfunction by h, then the theory in the Godel spacetime 
with a metric form 

d 2 s = (dt + e ax dy) 2 - dx 2 - \e 2ax dy 2 - dz 2 , (2.4) 

we have the equation 

h tt - 4e~ ax h ty + 2e~ 2ax h yy + h xx + ah x + h zz - M 2 h + rjh = 0. (2.5) 



Assuming the following form for the solution 

h(t, x) = expiikyij + ik z z — icut)^(x) (2.6) 

we can find from (2.5) find the differential equation of fy(x) 

+ a^, x - [u 2 + k 2 z + 2k 2 y e~ 2ax + Ak 2 y ue~ ax + M 2 — 77]* = 0. (2.7) 
If we introduce a new variable 

u = oT x lJ%k y \e- ax , (2.8) 
for the case of k 7^ 0, then (2.7) becomes 



l + ^ + ^M±^|* = , (2.9) 
1 4 au a l u l 1 



with e = k y /\k y \, which is the Whittaker's equation and its solutions are well known. 
From the asymptotic behavior of the solutions one can see that the condition to have an 
everywhere-bounded solution is 

r] = k 2 z + M 2 + a 2 [n +h + \a 2 - [u + V2ea(n + i)] 2 (2.10) 

with n =0, 1, 2, .... Note that the scalar wave equation in the Godel Universe has been 
solved in [21] and [22], where the quantized dispersion relation (Let i]=0 in (2.10)) was first 
derived. 
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When k y = 0, we can define a new variable 

v = e~ ax , (2.11) 

then (2.7) becomes 

2 7 2 

^«-^P* = 0, (2.12) 

The solutions are of simple powers of ,v and not bounded everywhere; thus we neglect it. 
To proceed, let us write 

Trln{A + M 2 ) = J d 4 xy/=g < t,x\ln(A + M 2 )\t,x > 

= f d 4 x^g~J2 < ^ X \ K >2< K \H A + M 2 )\ K >, (2-13) 



K 



where 



< t, x\K >= (2n)- 3/2 ^ L Q-2V2\k y \e- ax ^ exp(ik y y + ik z z - iut), (2.14) 
E = (-9)^ I "duo f dk y f dk z J2, (2.15) 

K n 

with K and L denote the sets of eigenvalues (u, k z , k y ,n) and (u, k z ,n) respectively. 

Two remarks are in order, (a) We need not worry about the existence of an integration 
over k y in (2.15) despite the fact that the eigenvalue is independent of k y . 

In fact, this integration over k y ,, could be used to perform the integration over the 
argument of ^ and thus makes the calculation of (2.13) possible, i.e. 

]T < t,x\K > 2 < K\ln(A + M 2 )\K > 

K 

= (-g)~ I duj dk z E < u; ' k z ,n\ln(A + M 2 )\u, k z ,n > 

J J n 

x(2tt)- 3 J dk y \^ L (^-2V2\ky\e- ax ^J \ 2 
= (2tt)~ 3 J did J dk z J2< uj,k z ,n\ln(A + M 2 )\iu,k z ,n> (2.16) 



5 



if the function \1> is normalized as / <i-u|\I>(-u)| 2 = I. (b) The weight {—g)^ 1 ^ 2 in (2.15) 
should be introduced because diudk is a tensor density rather than a tensor. Notice that 
we introduce a factor (— g)^ for dx, but {—g)~ for doodk z . This is because the weight of 
the density for a contravariant vector is the inverse of that for a covariant vector [23]. (The 
function exp(ik y + ik z — iui) in (2.6) tells us that (t, x) and (cu,k) are the contravariant 
vector and the covariant vector respectively.) It is the introduced factor (—g)~ which will 
cancel the factor e ax coming from the change of the variable from k y , to u, we can therefore 
have an x-independent effective potential. (Note that the Godel Universe is homogeneous in 
both time and space.) It seems that several previous articles have neglected this point. 
We therefore obtain 

Ta r 

V *ff = / ^EE < k z ,n,l\ln(A + M 2 )\k z ,n,l > (2.17) 

07T J n I 



where we have used the following replacement 

u + V2ea(n + -) -> i2irlT, (2.18) 

in order to perform a computation at finite temperature T. The expression (2.17) is divergent 
and needs to be regularized This is the work of section 3. 



3 The (-function regularization 

We now use the method of ( function regularization [20] to find the renormalized value of 
V e ff in (2.17). Introducing a complex parameter u, we obtain a relation 

/^EE < k z ,n,l\ln(A + M 2 )\k z ,n,l>= -—((v)\ u=0: (3.1) 
J ~ i dv 

where 

C(v)= /^EE<^V|(A + M 2 )|A; z ,n,/>^. (3.2) 

J n I 

Using (2.10) and (2.18) we have 

cm = / dk z e E^ 2 + m 2 + "V + \f + W + ^ 2 i 2 T 2 r 

n I ^4 
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r(i/2)i> 1/2) ^ ^ [m2 + ^ + i )2 + i_ a2 + 47r2/2Tr , 
r(i/2)r(i/-i/2),„/ - , i , 2 , i 2 , , o,^ 1 / 2 - 



r(iz) 



(E E ( M ' + -« 2 ^ + 7« 2 + 4vr 2 / 2 T 2 ) 
n j V 4 4 / 



-EE( m2 + « 2 ^ + 7« 2 + 4 ^ 2t2 ) 1/2 ")■ ( 3 - 3 ) 

n l V 4 / 

Although the above summation is formally divergent, it could be well defined through the 
analytic continuation with the formula [24, 25] 

J2(M 2 + (n/a) 2 + (l/b) 2 )^ = ab(M 2 y-»7rr(v) (l> - 1) + 2£'Z"- 1 tf 1 _ I/ (2Z) j , (3.4) 

n,l \ n,l J 

where K is the modified Bessel function and Z is defined by 

Z = TiM{n 2 a 2 + l 2 b 2 } 1 / 2 . (3.5) 

A prime on the summation denotes that the term n = I = is omitted. Using this formula 
(3.3) becomes 

((v) = (2naT)-^ ) n^(\a 2 + M 2 )f 2 ^ 



x 



_It> - 3/2) + £'[2Zr 3/2 K 3 /2-,(2Z) - Z»- 3/2 K 3/2 „ u (2Z) ] , (3.6) 

n,l 



where 

1/2 



Z x = 7T (i^a 2 + M 2 )(/ 2 /4tt 2 T 2 + 4n 2 /« 2 )) 
Z 2 = tt ((^a 2 + M 2 )(/ 2 /4vr 2 T 2 + n 2 /a 2 )) ^ (3.7) 



As T(z/) — > 1/z/ we obtain from (2.17) the final result 



V eff = (IGtt 3 )- 1 !^^ 2 + M 2 )f 2 ^-T,'[^r 3/2 K 3/2 ^Z) - Zr 3/2 K 3/2 ^2Z)]^j . 

(3.8) 



A comment should be mentioned before applying the above exact formula to analyze 
the cosmological phase transition. From the formula of analytic continuation (3.4) one 
sees that the cases of a = and/or T = will not make sense. However, (3.8) will give 
V e ff — > (l/247r)M 3 , as each term in the summation will approach zero if a and/or T — > 
(see below), which is in contrast to the known result. This difficulty can be solved by 
observing that because in the zeta regularization method the renormalization condition 
is not an a priori assumption, we may regard the result (3.8) as being 'renormalized' by 
V e ff = (l/24ir)M 3 at a and/or T = 0. Therefore one can obtain a useful result from (3.8) only 
if T n 1. This is because when T^I the temperature- independent terms, which are absent 
from (3.8), are much less than the temperature-dependent terms, which can be obtained 
from (3.8). We will separately discuss two cases: a « 1,T >> 1 and a » 1,T >> I. 

3.1 a « 1, while T » 1 

When a « 1 then Z 1: and Z 2 defined in (3.7) are very large. Using the properties 
K 3 / 2 (Z) — > ^JitZ~ 1 l 2 e~ z as Z » 1, we see that V e ff in (3.8) is dominated by the summation 
terms of n — ±1. Therefore we have the relation 

]r'[2Zr 3/2 K 3 /2-,(2Z) - Zr 3/2 K 3/2 ^(2Z)] 

n,l 

= -E' 4 + M 2 )(l 2 /4n 2 T 2 + l/a 2 )} 1 / 2 )^ 2 

i=i 

xK 3/2 (27r[(a/4 + M 2 )(/ 2 /4vr 2 T 2 + 1/« 2 )] 1/2 ]) 
= -(2/7r) 3/2 (a 2 /4 + M 2 y\-a 2 ex V \-2Tx{a 2 j ^ + M 2 )/a] 
+ ^(/ 2 /4vr 2 T 2 + l/a 2 )- 1 e^[-27r[(a 2 /4 + M 2 )(/ 2 /4tt 2 T 2 + l/a 2 )] 1/2] , 

1=0 

if a«l. (3.9) 

From this relation one sees that V e ff = when a = as mentioned in the above comment. 
When T >> 1 the summation over I in (3.9) could be approximated by an integration. After 
neglecting the T-independent terms the equation of (3.9) becomes 

-(2/7r) 3 / 2 (a 2 /4 + M 2 )- 1 J dl{l 2 /Aix 2 T 2 + l/a 2 )' 1 
xexp[-2n[(a 2 /4 + M 2 )a(l 2 /4n 2 T 2 + l/a 2 )] 1/2 ] 
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= -(2/7r) 3 / 2 (a 2 /4 + M^litaT [ dX—t—exp[-A{X 2 + lf% if a « 1, T » 1, 

J X + 1 

(3.10) 

A=^(^! + M 2 ) 1 / 2 . (3.11) 
a 4 

Substituting the value 

/•oo I 

/o rfx ^Ti e;rj9hA(x2 + 1)1/2] ^ (7rM)1/VA ' i f A » 1 > ( 3 - 12 ) 

which is evaluated in the appendix, into (3.10) and combining the tree level effective potential 
we obtain the final result 

V eff = m 2 + ^R- (2vr)- 2 a 3 / 2 T(^ + M 2 ) l ^exp{-^(^ + M 2 ) 1 ' 2 }. (3.13) 
Taking the limit <fi — > from the above equation we obtain the effective mass 

m eff = m 2 + ZR + A(2tt)-V/ 2 T + m 2 + 6R|" 1/4 - ^ + m 2 + 6R|~ 3/4 ) 

xexp ^-^|^ + m 2 + e J R| 1/2 ^ • (3.14) 

From the above equation we can find the critical temperature by setting m e //see that the 
symmetry will be restored at high temperature. 

3.2 a » 1 and T >> 1 

When a >> 1 and T >> 1 then Z\ and Z 2 defined in (3.7) are very small if n and 1 are not 
too large. In this case the summations over n and I in (3.8) could be approximated by the 
integrations. Also, because K 2 /s(Z) — > Z~ 3 / 2 — \Z 1 ^ we see that the main contributions 
in (3.8) are those coming from the summations over small n and 1. Therefore we have the 
relation 

£'[2Zr 3/2 K 3 /2-,(2Z) - Zr Z/2 K V2 _ v {2Z)] 

n,l 

9T / 2 \ -3/2 / 2 \ -V 2 

^ • M 2 S 1 + V2Tar- + M 2 ) S 2 (3.15) 



7T 2 V 4 / 1 V 4 
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where 



SAT, a) = [ dX [ dY(X 2 + F 2 )- 3 / 2 - / dX ( dY{X 2 + F 2 )- 3 / 2 

JXo J2Y JX JY 

S 2 (T, a) = [ X dX [ Y dY{X 2 + Y 2 Y 1 ' 2 - (* dX F dY{X 2 + Y 2 )' 1 ' 2 (3.16) 

JXo J2Y Jx JY 

X = 1/2T, Y = l/a (3.17) 

while X, and Y, being independent of T and a, stand for the suitable values corresponding 
to the reasonable cut-off values of n and I. Notice that the factor (a 2 /4 + M 2 )~ 3//2 in the first 
term of (3.15) will be canceled by the factor (a 2 /A + M 2 ) 3 / 2 in equation (3.8). Therefore 
the first term in (3.15) merely contributes a ^-independent term to the effective potential, 
which at most shifts a value to the effective potential, and we can drop it. Combining the 
tree level effective potential we finally obtain the result 

V eff = [m 2 + £R+ l\(27r)-^ 2 aTS 2 }<j ) 2 . (3.18) 
8 

From the above equation we also see that the symmetry will be restored at high temperature. 
Note that the critical temperature can be easily obtained ifT >> a » 1, because in this 
case we may set X = and thus S 2 is only a function of a. 

Note that in the case a 2 /4 + m 2 + £R — > then (3.14) seems to give a divergent m e ff, 
which is intuitively an impossible occurrence. The reason is that if a 2 /4 + m 2 + £R — > 
then Zi and Z 2 — > 0, (if <fi also approach zero), this corresponds to the case of a,T » 1. 
This means that for the models with several values of a, £ and m 2 (note that R = —a 2 ) 
which makes a 2 /4 + m 2 + £R — * we shall have the effective potential described by (3.18). 

4 Discussions 

After the new scenario of the inflationary Universe [2, 3] was proposed, many authors have 
been attracted to investigate the quantum-field effects on the cosmological phase transi- 
tion [7, 8]. The problem with the anisotropic effect on the Coleman- Weinberg symmetry 
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breaking has been studied [10-12]. We have also in a recent publication [ 13] studied the 
inhomogeneous effect on the Coleman- Weinberg symmetry breaking. 

Recently, Chimento et al [ 14] calculated the 1-loop effective potential for the scalar field 
in a Godel spacetime, and found that whether the symmetry was restored or not would 
depend only upon the curvature coupling constant £, irrespective of the magnitude of the 
rotation. 

As the early Universe is also characterized by non-zero temperature [17, 18] it is worth- 
while to study the temperature effect on the cosmological phase transition in a rotating 
spacetime. In this paper we have calculated the finite temperature 1-loop effective for a 4 
theory in the Godel spacetime. From our results one can find the critical temperature for 
the symmetry restoration, which in general depends on the curvature coupling constant and 
magnitude of the spacetime rotation. 

Finally we want to make a remark. As Godel spacetime is not very physical our cal- 
culation is not intimately related to the direct observation. However, the investigation in 
this paper will be (at least, we hope) helpful in the future analysis of the finite-temperature 
cosmological phase transition in the more general and more realistic (physical) spacetimes. 



Appendix 

Calculation of the integration in (3.12) 

roo , , roo 

J dX(X 2 + l)- 1 l 2 exp[-A{X 2 + l) 1 ' 2 } = J dY(Y 2 - xy^y^e'^ 

= e~ A / dZ(Z + \y l \Z(Z + 2)}- 1 ' 2 e- ZA 
Jo 

roo 

= e~ A dZ\Z{Z + 2)}- 1 ' 2 e~ ZA 
Jo 

roo 

-e- A / dZZ^iZ + lYHZiZ + 2)}- l l 2 e- ZA 
Jo 

The existence of an exponential factor e~ AZ in the integrand tells us that the dominant 
contribution of the integration comes from the neighbourhood around Z = when A » 1 
(i.e. a « 1, see the relation (3.11)). Because the first integral has a factor Z^ 1 / 2 while the 
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second integral has a factor Z 1 / 2 , the first integral must be larger than the first one. The 
first integral can be found exactly [26] and we have the result 

/ dX(X 2 + l)-^ 2 exp\-A{X 2 + If' 2 } = K (A) 
Jo 

where K is the modified Bessel function. Taking the expansion of K for A >> 1 the equation 
(3.12) is obtained. 
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